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The  work  accomplished  under  this  contract  was  directed  toward  the 
development  of  a theoretical  computer-oriented  mode]  of  the  polar 
ionosphere  between  the  heights  of  50  and  500  km.  This  model  was  based 
on  first  principles  of  atmospheric  physics,  except  where  the  physical 
processes  were  not  understood.  The  model  was  to  be  comprehensive 
enough  to  predict  both  the  "quiet"  and  "disturbed"  polar  ionospheric 
conditions.  Our  strategy  for  solving  these  problems  was  to  model  the 
"quiet"  ionosphere  first  by  separating  the  total  height  region  into 
the  E and  F regions  (100-500  km)  and  the  D region  (50-100  km).  This 
was  done  because  the  physical  processes  are  much  different  above  and 
below  ^00  km. 

E and  F Regions 

The  E and  F regions  of  the  polar  ionosphere  are  similar  to  those 
at  raid-latitudes  except  for  energetic  particle  precipitation  down  the 
geomagnetic  field  lines.  To  handle  the  kinetic  processes,  we  modified 
a chemical-kinetics  computer  program  developed  by  Adams  and  Megill 
[1970].  This  computer  code  allows  as  many  as  31  species  to  interact 
in  up  to  51  different  reactions  with  15  separate  photorr act  ions . The 
flow  diagram  of  the  modified  computer  program  is  given  in  Figure  1. 

The  processes  indicated  in  the  dashed  areas  were  not  implemented  before 
termination  of  the  contract. 

Photoionization  rates  during  the  daytime  for  0+,  N2+,  0+,  and 

+ 

NO  are  calculated  by  utilizing  the  approximations  to  the  Chapman 
function  for  optical  depth  derived  by  Swider  [1964],  During  the  daytime 
the  photoionization  rates  are  three  orders  of  magnitude  greater  than 
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the  ionization  produced  by  the  precipitation  "drizzle"  electron  flux. 

See  Figure  2.  The  nighttime  ionization  in  the  polar  ionosphere  is  a 

complicated  mixture  of  energetic  particles  and  scattered  solar  radiation 

/ 

The  mid-latitude  E region  is  maintained  at  a level  of  10-10*  electrons/ 

3 

cm  through  the  night  by  solar  radiation  scattered  to  the  dark  side  of 
the  earth  by  the  earth's  hydrogen  and  helium  geocoronia.  The  model  of 

Ogawa  and  Tohmntsu  [1967]  was  utilized  even  though  the  ionization  rate 

from  scattered  light  is  a function  of  the  radiation-transport  geometry 

[Meier,  1974].  Comparison  of  the  production  rates  due  to  scattered 

radiation.  Figure  3,  and  drizzle  flux.  Figure  2,  show  that  these  two 

processes  are  about  the  same  in  the  nighttime  polar  ionosphere.  F.lectron 

precipitation  was  divided  into  a quiet  "drizzle"  flux  and  into  3 

strengths  of  auroral  arc  precipitation.  The  "drizzle"  precipitation 

is  that  found  in  the  quiet  auroral  oval  and  in  the  polar  cap  [Burch, 

1968].  The  auroral  arc  precipitation  was  designated  TBC-I,  II  and  III, 
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and  defined  as  having  flux  values  of  10  exp(-  K/3.5,  10  exp(-  E/8.5), 

and  10  exp (-e/12.5)  electrons  per  cm^-sec-stcr-kev  where  F is  the 

particle  energy  in  kilo-electron  volts.  The  production  rates  for  N^+, 

■f  *4" 

0^  and  0 due  to  electron  fluxes  is  given  by 


Qi  (h) 


°iNi 


o N + o J + a N • 
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Q(h)  (cm-3  sec-1)  i=N  ,0  or  0 

^ 2 


where  Q(h)  is  the  Integra1  over  all  energies  of  the  ionization  rate 


Production  Rates  from  Solar  Radiation  and  "Drizzle"  Precipitation  Electrons 
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q(!i,F)  given  by  Lazarev  1 1 967 ) at  the  penetration  depth  X (K)  deduced 
by  Gerard  [ 1970 J 

, to 

= j q(h,K)  dl-  (cm  ^ sec  ^). 

*'min 

The  Ionization  rate  for  an  IBC  II  auroral  arc  is  illustrated  in  Figure 
3. 

Utilizing,  the  ionization  sources,  a reaction  set  such  as  given 
by  Rees,  et  al,  [1967],  and  the  neutral  atmosphere  model  of  Jacchia 
[Jacchia,  1971],  it  is  possible  to  calculate  the  time  dependent  densities 
of  the  ionic  species.  The  computer  code  was  very  expensive  to  run  on 
the  Utah  State  University  Burrough's  Computers,  therefore,  only  a 
limited  amount  of  computation  was  actually  done.  Figure  4 illustrates  the 
nighttime  ionic  densities  at  2300  local  time,  65 . IN . latitude ,212. 5°W  longitude 

produced  by  scattered  radiation  and  drizzle  flux.  Also  shown  are  ion  densities 
after  20  minutes  of  an  auroral  described  by  the  flux-energy  diagram 
of  Figure  5.  The  0+  density  is  excessively  high  at  300  km  because 
the  calculations  neglected  vertical  transport  by  diffusion. 

The  production  rates  for  N2+,  0+,  and  0+  permit  the  calculation 
of  the  volumn  emission  rates  for  the  Nitrogen  3914>\,  the  oxygen  5577A 

O 

and  6300A  radiations  (Rees  et  al,  196-7] . The  3914A  emission  results 
from  electron  impact  ionization  and  excitation  of  N^.  This  emission 
is  imoortant  because  of  the  proportionality  between  the  primary 
electron  flux  and  the  N2  production  rate.  The  oxygen  green  and  red 


Theoretical  ion  densities  for  quiet  (Q)  and  auroral  arc  (A)  conditions 
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linos  are  produced  by  the  reactions. 

0(!S)  -*•  0(]D)  + hv  (5677A) 

0( 1 D)  » 0(31‘)  + hv  (6300A) 

The  0(  S)  state  is  populated  by  collisionnl  excitation  of  secondary 
electrons  with  oxygon  atoms  and  also  the  dissociative  recombination 
of  . The  0(^D)  state  is  populated  by  collisions  with  the  ambient 
electrons  as  well  as  the  same  excitation  processes  for  the  0(^S)  state. 

Our  calculation  of  emission  rates  is  based  on  the  model  of  Rees  ct  al , 
[1967],  and  includes  collisional  quenching  but  not  diffusion.  Figure 
6 is  a plot  of  the  three  radiation  profiles  for  quiet  and  disturbed 
(the  flux  of  Figure  5)  conditions.  This  is  as  far  as  the  work  in 
this  area  progressed  before  the  contract  was  terminated.  Additional 
effort,  funded  by  the  University,  is  being  made  to  evaluate  this 
computer  code  by  comparing  simultaneous  electron  fluxes  and  radiation 
profiles  which  were  measured  by  a research  group  at  Utah  State  University. 
Copies  of  the  completed  work  will  be  sent  to  the  distribution  list 
at  the  time  of  completion. 

Many  of  the  problems  encountered  in  developing  a polar  ionospheric 
model  were  at  the  forefront  of  current  research.  Such  was  the  case 
in  dealing  with  motions  produced  by  the  collisional-heating  of  precip- 
itating particles  and  auroral  electric  fields,  and  in  simply  but 
adequately  describing  the  chemistry  of  the  D region  (below  100km). 

The  Joule  heating  produced  by  auroral  electric  fields  is  of  the  same 


photons c®  £cc 

Figure  6.  Volumn  emission  profiles  for  quiet  (Q)  and  auroral  arc  (A)  conditions 
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order  of  magnitude  as  the  solar  ultraviolet  heating  of  the  thermosphere 
and  results  in  an  important  energy  sourro  in  the  polar  regions  especially 
during  magnetically  disturbed  conditions.  Since  the  air  is  easily  set  into 
motion,  we  would  expect  this  heating  to  cause  large  winds  in  the  polar 
ionosphere.  The  problem  of  describing  the  auroral  motion  was  begun 
by  modeling  a simplified  2-dimensional  nor Lh-south-ver tical  slab  of  the 
auroral  zone.  Motion  was  caused  by  a north-south  electric  field  driving 
ions  across  the  geomagnetic  field  lines  and  imparting  heat  energy  and 
momentum  to  the  neutral  particles  through  collisions.  The  gas  equations 
for  mass  continuity,  momentum  and  energy  including  viscosity  and  heat 
conduction  were  linearized  by  assuming  small  amplitude  (perturbation) 
motions.  A vertical  pressure  coordinate  was  employed  together  with 
the  assumption  of  hydrostatic  equilibrium.  Motion  in  the  east-west 
direction  was  calculated,  but  no  east-west  variations  or  pressure 
gradients  were  allowed.  The  neutral  gas  motion  produced  by  a 50  millivolt/ 
meter  electrical  field  acting  on  a static  resting  atmosphere  for  3-1/2  hours 
is  illustrated  in  Figure  7.  This  important  work  is  not  completed,  but 
Dr.  Harris  has  received  free  computer  time  at  the  National  Center  for 
Atmos,  heric  Research  to  make  additional  calculations,  especially  of 
motions  about  an  arc.  If  desired,  a preprint  of  this  paper  will  be 
sent  to  the  distribution  list  when  available. 

A complete  chemical  model  of  the  D region  is  still  under  development 
due  to  the  multitude  of  possible  chemical  reactions.  Our  philosophy 
for  modeling  this  region  was  to  utilize  equivalent  reaction  rates  that 


Fleur*  7.  Wind  profiles  generated  by  auroral  rone  electric  field 
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were  compatible  witli  observed  ionization  features  so  that  a predictive 
model  could  be  implemented  without  having  to  evaluate  the  large  reaction 
set  that  is  in  current  vogue.  The  first  problem  encountered  with  this 
model  was  the  calculation  of  the  solar  ionization  and  dissociation 
rates.  We  found  that  tiie  NO  ionization  rate  due  to  the  1027-1337  A 
solar  flux  region  had  not  been  calculated  utilizing  detailed  structure 
of  either  the  02  absorption  cross-sections  or  'he  NO  ionizatio.  cross- 
seclions.  A computer  subroutine,  incorporating  the  details  of  the 
measured  N'0  ionization  cross  sections  [Watanabe  et  al , 1967],  and  the 
02  absorplion  cross  sections  [Adams,  1974],  was  written  to  evaluate 
ionization  in  Lhis  hand.  Results  for  an  overhead  sun  are  shown  in  Figure 
8.  It  is  seen  that  the  hydrogen  Lyman  a line  (1215. 7A)  produces  almost 
an  order  of  magnitude  more  ions  than  the  1027-1 337A  solar  band.  Thus, 
this  band  is  only  of  minor  importance  in  the  D region  chemistry. 

The  dissociation  rates  per  molecule  for  Cll^ , CO,  NO,  C02 , H20, 

N20,  NO^,  ^2H^,  H202,  02  and  0^  have  been  calculated  versus  altitude 
and  zenith  angle  for  solar  wavelengths  between  1000A  and  the  dissociation 
threshhold.  Absorption  cross  sections  for  these  species  were  obtained 
from  the  Joint  Institute  for  Laboratory  Astrophysics  (reference  library) 
of  the  University  of  Colorado,  Boulder,  Colorado.  The  absorbing 
atmosphere  was  assumed  to  be  oxygen  (Jacchia,  71  and  U.S.  Standard 
Atmosphere,  62)  and  Ozone  [Adams  and  Megill,  1970],  These  calculations 
arc  being  prepared  for  publication  in  a technical  journal  and  represent 
the  most  comprehensive  and  detailed  description  of  dissociation  rates 
available  today.  Figure  9,  from  this  paper,  shows  the  altitude  where 


